Evolutionary Applications

1Effective population size (Ne) is a critical parameter for evaluating the evolutionary and persistence potential of endangered populations and for designing sustainable conservation strategies. Captive breeding and release programs are widely used across taxa to reduce risk of extinction when natural reproduction is insufficient or no longer possible, making it essential to assess their consequences. We used the case study of the landlocked Saimaa salmon (Salmo salar), one of the most critically en-dangered salmonid populations in Europe, with unique evolutionary significance due to its isolation from other populations since the last glaciation. Using long-term demographic data (1969-2024) from wild-caught founders of a captive breeding and release program, we estimated the effective population size under multiple scenarios of variance in reproductive success. Across scenarios, Ne ranged from 33 to 81 individuals, representing 32%-75% of the census size. Captive breeding practices aimed at equalizing parental contributions during fertilization and early life stages increased Ne by 12% compared to natural reproductive conditions. However, variation in survival after early developmental stages, typically beyond direct management control, remained a key determinant of Ne. Despite recent increases in the number of founders, the population remains genetically vulnerable due to historical bottlenecks. These results highlight that while captive breeding programs can partially mitigate genetic risks, their effectiveness depends critically on both controlled and uncontrolled sources of variance in reproductive success. Strengthening such programs may require combining breeding management with habitat restoration and, where appropriate, genetic rescue to ensure the long-term evolutionary potential of such unique and endangered populations.

Genomic erosion as a manifestation of small effective population size (Ne) and consanguinity subverts long-term perpetuation of threatened species by compromising their adaptive potential; however, the integration of genomics remains limited in applied conservation efforts to guide priorities. This study combines non-invasive sampling, double-digest Restriction site-associated DNA sequencing (ddRAD), and population-genomic analyses to assess genetic health in two vulture assemblages-mixed wild enclosure and captive breeding cohorts. Both the geographical locations exhibit signs of populations in distress: low genetic diversity and abundant intermediate-length runs of homozygosity (RoH), consistent with long-term reduced Ne plus recent demographic isolation. Our demographic model runs favoured ancient migration (AM) topology characterised by an ephemeral window of gene flow, taken over by a prolonged population separation period. The mutation quantification results from approximately 59,000 outgroup-polarised SNPs reveal higher additive burden and more homozygous-derived sites in BKN. However, this was later traced to low-impact and non-coding variants rather than a surge in the loss-of-function (LoF) alleles. The data support a genomic profile that carries an elevated risk from polygenic/aggregate deleterious burden in BKN despite a scarcity of high-impact mutations. By highlighting the disconnect between genetic resilience and demographic recovery, our results accentuate the need to incorporate genomics-informed inbreeding and monitoring programs, while also focusing on reducing anthropogenic mortality with genetic augmentation.

Genomic forecasting approaches based on genotype-environment associations (GEAs) are increasingly used to estimate genomic offsets (GOs), which predict population maladaptation and extinction risk under current or future climatic conditions. Despite their widespread use, only a subset of studies have evaluated how accurately GOs predict (mal)adaptation, limiting their interpretation and application in policy and management. Here, we used GEA analyses to estimate GOs for past, present, and future climates in Lycaeides butterflies, focusing on the causes of variation in GOs among populations and their relationships with demographic parameters inferred from population genomic data. Using multivariate linear regression and genotyping-by-sequencing data from 42 Lycaeides populations (922 butterflies), we found that mean annual temperature, cumulative annual precipitation, and hybridization history together explained 47.6% of variation in genome-wide allele frequencies. Genomic offsets differed substantially among populations and across past, present, and future climates, with evidence for increasing maladaptation under more distant future climate scenarios. We found no relationship between GOs for present climates and contemporary effective population size. In contrast, genetic diversity, which reflects long-term effective population size, and local rates of gene flow together explained 27.3% of variation in contemporary GOs. Populations with higher genetic diversity and more gene flow exhibited lower GOs, consistent with the hypothesis that genetic diversity enhances adaptive capacity and that gene flow may introduce adaptive alleles. Overall, our results support the utility of GO predictions, particularly when validated with independent measures of adaptation, while cautioning against simplistic interpretations of GO as a direct measure of maladaptation in conservation and management contexts.

The fitness of immigrants and their descendants determines the effectiveness of gene flow. Genetic incompatibilities or outbreeding depression can limit the spread of novel alleles, while highly fit immigrant lineages can hasten introgression. These fitness effects of gene flow can also differ between generations as immigrant and resident haplotypes recombine. Understanding the genetic factors that shape immigrant fitness over multiple generations is increasingly important as habitat fragmentation threatens populations by reducing genetic variation and leading to increased levels of inbreeding. Few studies have measured the multigenerational fitness of immigrant lineages, especially within populations that had histories of high gene flow. We used 33 years of life history and pedigree data on a population of Florida scrub-jays (Aphelocoma coerulescens) with historically high immigration to quantify the fitness of immigrants and their descendants. We compared the fitness of immigrants and residents as well as their resulting descendants (F1, F2, etc.) to determine the composite genetic effects responsible for fitness differences. We found evidence of additive benefits of immigrant ancestry and heterosis driven by non-additive effects that persists for multiple generations. These results are promising for conservation efforts aiming to increase connectivity and illustrate the complex dynamics that determine the rates of introgression in natural populations.

Allopatric isolation under contrasting environments can drive rapid phenotypic divergence, even over contemporary timescales. Rapid changes in morphology or physiology can allow organisms to adapt to biotic and abiotic characteristics of their habitats. While studying metabolism, growth and resources needs may allow to understand adaptation to several selective pressures, these traits are rarely jointly considered. We investigated morphological, growth, and metabolic divergence in two allopatric populations of Arctic charr (Salvelinus alpinus) sharing a common evolutionary origin but inhabiting contrasting environments. We combined field observations, common garden and quantitative genetic approaches to disentangle contributions of genetic divergence and plasticity to phenotypic variability. Wild adults differed in body shape and growth trajectories, potentially reflecting plasticity related to resource availability and temperature variations. Under common garden conditions, juveniles displayed inter-population differences in routine metabolic rate, its allometric scaling with body mass. These patterns suggest divergent selection on physiological traits. Despite low neutral genetic differentiation, phenotypic divergence unfolded in fewer than 100 years, suggesting that plasticity and selection can promote rapid multi-trait changes. These findings highlight that considering changes in physiological, growth and morphological traits can reveal the adaptive potential of small, isolated populations facing rapid environmental change.

RDL (Resistance to dieldrin) is a GABA-gated chloride channel that was first described as target of the insecticide dieldrin. Despite dieldrin being discontinued for decades because of its environmental per-sistence and health concerns, Rdl resistance mutations (A296S, A296G) continue at high frequencies in natural populations of the malaria mosquito Anopheles gambiae complex across Africa, suggesting a selective advantage. We have recently shown that RDL acts as a critical modulator of mosquito auditory sensitivity. Because acoustic recognition is essential for mate acquisition in An. gambiae, we hypothesized that these mutations confer a pleiotropic effect on mating success in the field, mediated through altered acoustic sensitivity, with potential consequences for sexual selection. We first provide laboratory evidence that resistance mutations enhance auditory behaviours of An. gambiae and show that the effect of environmental noise on mating success depends on the male Rdl genotype. We then conduct field collections in the city of Bangui (Central African Republic) and surrounding rural areas, revealing the presence of Rdl resistant alleles and their association with the urban environment, and within the city, with the noisiest locations. We also show decreased mating success of susceptible females with increasing noise levels, suggesting detrimental effects. Together, our findings support that Rdl resistance mutations enhance auditory function and mating success in acoustically challenging environments. We propose that this auditory advantage may contribute, together with other selective pressures such as cross-selection by other insecticides, to the persistence of these alleles in nature and may facilitate urban colonization by malaria vectors. Our study reveals, for the first time, an unintended evolutionary consequence of insecticide use, where a resistance mutation has been co-opted to enhance sensory performance and ecological adaptation, with significant implications for vector management strategies.

O_LIEffective conservation of highly migratory species requires understanding genetic structure across breeding populations and access high{square}resolution markers capable of assigning individuals from mixed aggregates (e.g. bycatch or new nesting sites) to their natal origins. Genomic approaches provide unprecedented resolution but add methodological challenges; thus, it is essential to first build a genomic baseline from known breeding areas and then evaluate strategies for assigning unknown individuals. C_LIO_LITo address this, we used 2b-RAD sequencing, a genomic reduction technique useful for degraded DNA, and loggerhead turtles as a case study. This species shows philopatric breeding, while juveniles and adults form mixed aggregations in foraging grounds. C_LIO_LIOur results highlight the importance of building baselines that include all potential source populations contributing to mixed aggregations. We detected hierarchical genetic differentiation and high resolution and successfully assigned the natal origin of 124 unknown individuals from four Mediterranean foraging grounds. These grounds showed distinct source contributions, and comparisons with previous studies suggest possible temporal shifts in stock composition. C_LIO_LIWe provide a comprehensive genomic baseline for individual assignment of Altanto-Mediterranean loggerhead turtles of unknown natal origin and a general framework for identifying population-specific threats in highly migratory species. C_LI

The sunflower sea star (Pycnopodia helianthoides) suffered a catastrophic population decline across its range from 2013 to 2017 due to the devastating Vibrio pectenicida FHCF-3 driven sea star wasting disease (SSWD) pandemic with minimal signs of population recovery. The functional extinction of this apex predator across substantial parts of its range has created a need to identify and track the remaining intact populations. Environmental DNA (eDNA) approaches provide a simple, cost-effective, and non-destructive method for monitoring occurrences, and in some cases abundances, of marine species, consistently outperforming visual occurrence monitoring efforts in sensitivity, speed, and cost. Here, we designed, developed, and validated a P. helianthoides-specific eDNA assay to identify refugia, using both quantitative and digital droplet PCR approaches. We first generated the most comprehensive sea star mitochondrial genome reference database to date (n=93 taxa, n= 15 novel). We then used unikseq and Geneious bioinformatics software to identify the unique nad5 gene region and design a highly specific hydrolysis probe-based PCR assay. We validated the performance of this assay through laboratory, mesocosm, and field testing, demonstrating a highly specific and sensitive assay. In a field application of the new assay across regions in British Columbia, Canada, we found a positive correlation between P. helianthoides eDNA concentrations and biomass density, especially when appropriately accounting for spatiotemporal integration scales (R2=0.67). The eDNA assay provides a rapid and scalable tool for monitoring the sunflower sea star which has been proposed for listing as threatened under the U.S. Endangered Species Act of 1973. Molecular tools like the one presented here enhance management and recovery efforts not only by identification and monitoring of remnant wild populations, but also by helping to assess population level response and recovery following reintroduction efforts.

The deployment of clothianidin-based insecticide formulations in malaria vector control has highlighted the capacity of Anopheles funestus to displace more susceptible mosquito species in treated areas and to rapidly evolve resistance under selection pressure. Metabolic detoxification, together with structural and genetic changes in nicotinic acetylcholine receptors (nAChRs), the primary molecular targets of neonicotinoids, can reduce insecticide efficacy. Here, we characterized amino acid substitutions across all 11 nAChR subunits in An. funestus to assess standing variation that may facilitate adaptive responses to chemical exposure. Using whole-genome sequencing data from 656 mosquitoes sampled in 13 African countries, we found marked contrasts in the distribution of nonsynonymous variants among nAChR subunits. Most subunits are strongly constrained and carry no missense variants, whereas two loci (3 and 7) display three geographically widespread amino acid substitutions across the continent. In contrast, 9 and {beta}2 accumulate dozens of nonsynonymous mutations occurring at intermediate to high frequencies, including within domains involved in orthosteric ligand binding and channel gating. Genetic differentiation at nAChR loci among populations from different countries is low to moderate, although several nonsynonymous mutations display high FST values consistent with geographic structuring. These results highlight relaxed constraint on two subunits that may provide opportunities for evolutionary diversification within a conserved family of multimeric receptor assemblies. Such diversification has not been observed in vector species displaced by An. funestus in indoor residual spraying areas, and the potential implications for reduced sensitivity to neonicotinoids are discussed.

Strong population contractions can leave a persistent genomic legacy that can influence populations long after their demographic recovery. While bottlenecks facilitate the removal of strongly deleterious mutations, the effectiveness of purging may be limited in historically small populations. The Kirtlands warbler (Setophaga kirtlandii) is a rare North American songbird with an ancestrally small population. After narrowly evading extinction, they are one of few species that have been delisted from federal protections in the USA. Despite their recovery, a previous study showed evidence for recent inbreeding and a high burden of deleterious mutations that may have not been purged despite strong bottlenecks. Historical DNA offers a unique opportunity to understand the consequences of recent demographic declines on genetic diversity. Here, we use DNA from over 100-year-old museum specimens to estimate changes in genetic load in the Kirtlands warblers pre- and post-bottleneck. We validate our results with forward-in-time genetic simulations and explore how sample size and missing data can affect estimates. Both empirical data and simulations suggest a reduced ability to purge deleterious mutations in this historically small population. Our simulations also highlight that limited sampling design and data quality can constrain the ability to detect changes.

Matching habitat choice provides a mechanism for individuals to maximise their expected fitness by selecting an environment that better fits their phenotype. Many animals choose their local environment by evaluating levels of perceived predation risk against possible resource gain. To test if predation risk is a major driver of habitat choice, we quantify scototaxis, or preference for dark versus light backgrounds, in juvenile guppies. As light backgrounds increase visibility to predators, this aspect of habitat choice captures variation in boldness in small fishes. By rearing and testing 586 fish descended from ten natural populations from Trinidad under common garden conditions, we first quantify (broad sense) heritable variation, i.e. evolutionary potential, within populations. Next, we test for evolutionary divergence among populations in mean preference, and if present, whether ancestral predation regime is a mediator of divergence. Finally, we ask whether families and/or populations differ in the amount of behavioural variation they contain. Habitat choice varied among families (12% of total variance), consistent with heritable variation (0.2). We also found mean preference varies among populations (11% of total variance explained). Evolutionary divergence among-populations is partly explained by ancestral predation regime, with populations from low-predation sites showing a stronger average preference for dark backgrounds than high-predation populations from the same river. Additionally, we find that within-population behavioural variation is greater in high-predation populations. We conclude that guppy populations contain heritable variation that could facilitate adaptive evolution if scototaxis is subject to natural selection. Furthermore, while genetic drift may also contribute to evolutionary divergence among-populations, observed patterns are qualitatively consistent with local adaption to predation regime. Our results suggests that high predation sites favour bolder habitat choice on average, but also that local predation regime shape the evolutionary dynamics of variation, perhaps by maintaining shy-bold variation among-individuals or by favouring individuals with less-predicable behaviour.

Telomere length (TL) is increasingly used in ecology as a biomarker of individual quality and environmental stress, yet research on non-model species with complex life histories remains limited. Because TL varies among tissues and across ages in a species-specific manner, identifying non-lethal tissues that reliably reflect whole-organism telomere dynamics is essential for longitudinal telomere studies in the field. This study aimed to evaluate tissue-specific TL in Japanese eel (Anguilla japonica), an endangered catadromous fish. We first mapped the chromosomal distribution of telomeric sequences using fluorescent in situ hybridization (FISH), the first application of this method in this species. We then tested whether muscle and caudal fin, which can be sampled easily and non-lethally, can serve as suitable proxy tissues for TL measurements in wild individuals. Relative telomere length (RTL) was quantified by qPCR in blood, brain, caudal fin, gonads, heart, liver, and muscle. FISH analysis confirmed telomeric repeats at all chromosomal ends, with only weak interstitial signals on three chromosomal pairs unlikely to affect qPCR-based estimates. A generalized additive mixed model and Wilcoxons signed-rank tests revealed significant inter-tissue differences: RTL was shortest in the brain and muscle and longest in liver, blood and caudal fin. Muscle and caudal fin RTL were significantly correlated with RTL in many other tissues, supporting their use as proxy tissues for longitudinal TL monitoring, including responses to environmental variation. Both total length and age were tested as explanatory variables for RTL, and the model including total length showed a better fit than the age-based model. Non-linear relationships between RTL and total length observed in several tissues suggest physiological shifts associated with growth and sexual differentiation. Overall, these findings advance understanding of telomere dynamics in eels and establish muscle and caudal fin as suitable tissues for repeated, non-lethal TL assessment in ecological and conservation contexts.

AO_SCPLOWBSTRACTC_SCPLOWEnvironmental heterogeneity across freshwater systems often promotes phenotypic variation, yet disentangling environmentally induced variation from heritable differentiation remains a central goal in evolutionary ecology. We investigated the geographic distribution and morphological differentiation, and heritability of shell traits among populations of the freshwater lymnaeid snail Pectinidens diaphanus in Patagonia. Extensive field surveys across 196 freshwater sites revealed that the species occupies a broad range of lentic and lotic habitats and constitutes the only lymnaeid inhabiting southern Patagonia. While reproductive anatomical structures were conserved across populations, shell shape differed markedly among populations from contrasting habitat types, with population identity explaining nearly 50% of total shape variation. Populations from hydrologically unstable habitats (ponds and streams) exhibited more elongated shells and relatively smaller apertures, a pattern consistent with functional responses to hydroperiod variability and desiccation risk. To assess the heritability of this differentiation, we conducted a common-garden experiment across two generations. Shell shape differences between permanent- (lagoon) and temporary- (pond) habitat-derived populations persisted into the G2 generation reared under standardized laboratory conditions, indicating that the observed variation is not solely a response to local environmental conditions but includes a heritable component. Together, our findings demonstrate that P. diaphanus constitutes the sole lymnaeid across southern Patagonia, occupying a broader range than previously documented, and that populations show heritable shell differentiation potentially associated with contrasting freshwater habitats. By integrating large-scale biogeographic surveys with morphometric and experimental approaches, this study provides new insight into how habitat variation may contribute to ecological and evolutionary differentiation in freshwater gastropods.

The Red Queen Hypothesis proposes that genetic variation is maintained in populations through antagonistic coevolution of hosts and parasites. A major assumption of the Red Queen Hypothesis is tight genetic specificity for infection. However, it has been argued that this genetic interaction of host and parasite (GHxGP) is sensitive to environmental context (GHxGPxE). Environmental change could accordingly disrupt coevolutionary oscillations on relevant time scales, calling into question antagonistic coevolution as a general and robust explanation for the maintenance of genetic diversity. To evaluate this critique, we used the plant-parasitic nematode Meloidogyne arenaria and its natural bacterial parasite Pasteuria penetrans to determine if specificity is altered by temperature. We exposed six isofemale host lines to five parasite sources at three ecologically relevant temperatures. We found that, at two of three temperatures, susceptibility to infection depended on the specific combination of host line and parasite source (GHxGP). This specificity varied across temperatures, consistent with a GHxGPxE effect. This three-way interaction was driven both by quantitative changes in the strength of specificity across temperatures and shifts in the susceptibility rankings of host-parasite combinations. Our study contributes a rare experimental test of a proposed challenge to the Red Queen Hypothesis and suggests the potential for environmental context to change host-parasite specificity.

Conspecific sperm precedence via cryptic female choice is a post-ejaculatory selection process that reduces hybridization, and can be pronounced in sympatric species. In their native Europe, Atlantic salmon (Salmo salar) and brown trout (Salmo trutta) exert conspecific sperm precedence under heterospecific sperm competition, which is at least partially enabled by female reproductive fluid. We examined post-ejaculatory selection of both species in Newfoundland, Canada, where Atlantic salmon evolved in absence of brown trout, but now experience hybridization threats due to anthropogenic introductions. Using split-ejaculate and split-clutch in-vitro fertilizations we evaluated whether allopatric evolution has relaxed this selection in Atlantic salmon, and found that they had no ability to bias paternity towards conspecific males, whereas naturalized brown trout retained a strong ability to do so. Female reproductive fluid influenced this, as when fluid associated with a species eggs was swapped, hybridization increased. In the artificial situation of no female reproductive fluid during sperm competition, paternity changed dramatically, but sperm swimming performance did not predict it. Our findings contribute to understanding the evolution of cryptic female choice and how the mechanisms of reproductive isolation can be reinforced through sympatry, while also highlighting a new potential conservation concern for North American Atlantic salmon.

Hybridisation between wild and domestic taxa can favour the spread of domestic alleles into wild populations through backcrossing. The complex interplay of random genetic drift, recombination, and selection can shape the fate of introgressed alleles. Maladaptive domestic variants are likely to be purged by natural selection, but others may persist across generations. It has long been known that the Apennine Italian wolf population, exposed to large numbers of free-ranging dogs, has experienced extensive introgression. The unusually high frequency of black wolves observed in Italy, compared to other European populations, may parallel patterns documented in North American wolves, where the melanistic KB allele at the CBD103 gene, of domestic origin, has spread over thousands of years of introgression. We tested whether the KB mutation entered the peninsular Italian wolf population via hybridisation and spread through adaptive introgression. Genome-wide analyses of black and wild-type (grey-coated) Apennine wolves showed no clear signatures of recent dog ancestry in most melanistic animals. Our ancestry reconstruction approaches identified two distinct KB haplogroups of domestic origin, suggesting multiple introgression events. Notably, we found molecular evidence consistent with balancing selection on the KB haplotypes, whose functional role, nonetheless, warrants further research. Therefore, the microevolutionary genomic and ecological consequences of wolf-dog hybridisation in Italy should be carefully investigated to inform appropriate science-based conservation management strategies.

In many tropical areas, seasonal rainfall leads to distinct dry and wet seasons. Many butterflies developing under wet season conditions develop into adults with large ventral eyespots on the wing margins, whereas those developing under dry season conditions have smaller or no eyespots. In greener, wet season habitats, larger eyespots can divert predator attacks toward the wing margins, while reduced eyespot size improves camouflage in the dry leaf litter-dominated habitat during the dry season. However, the dry season is also characterised by higher desiccation stress than the wet season. We hypothesised that larvae developing under dry season conditions develop into adults with higher desiccation tolerance than those reared under wet season conditions. We tested this by rearing larvae of the butterfly Mycalesis mineus under simulated dry and wet season conditions and assaying the desiccation tolerance of the resulting adults. Butterflies reared in dry conditions survived longer under desiccation stress, lost lesser water during pupal-adult metamorphosis, and were heavier than those reared in wet conditions. We also tested the correlation between eyespot size and desiccation tolerance. A negative correlation between the traits would be expected if similar developmental pathways regulate them. Consistent with this expectation, individuals with smaller eyespots had higher desiccation tolerance. Our results demonstrate plasticity in desiccation tolerance, and suggest that predator avoidance and desiccation tolerance traits may share similar developmental pathways.

Climatic and biogeographic variables are often used as a proxy for tree genetic diversity, but local factors can also influence it. We propose that woodland age, presence of ancient trees, and population size could impact genetic diversity. Using the RBG Kew UK National Tree Seed Project as a study case, we examined how these factors are accounted for during seed collection. We found 42% of tree seed collections come from ancient woodlands and that 8.4% overlap with ancient trees. Sampled forest patches size ranges from few individuals to several thousand. We then carried out a pilot to examine the role of population size on functional traits variation, testing the relationship between population size and seed germination and seedling thermal stress sensitivity in three populations of the Betula pubescens Ehrh. complex. We found that seeds and seedlings from larger populations showed higher fitness and stress resistance. Our results highlight the importance of local factors to predict variation in functional traits, relevant for tree resilience. Existing seed collections of native species stored in conservation seed banks offer a valuable resource to explore these factors and improve our understanding of genetic diversity in tree populations, with implications for biodiversity conservation and forestry production.

The marbled teal (Marmaronetta angustirostris) is a widely distributed but declining waterfowl species, classified as Near Threatened globally and Critically Endangered in Spain. Despite ongoing conservation actions, including ex situ management and population reinforcement programmes, the genomic consequences of long-term captivity, inbreeding, and patterns of functional genetic variation remain unknown due to the absence of a species-specific reference genome. Here, we present the first chromosome-level genome assembly for this species. The genome was generated using PacBio HiFi long reads and Omni-C data, yielding a 1.15Gb assembly with a scaffold N50 of 76.95Mb. A total of 97.16% of the assembly was anchored into 36 chromosome-scale scaffolds, including the Z and W sex chromosomes. BUSCO analysis recovered 99.2% of conserved avian genes. Gene prediction was performed using both ab initio and homology-based strategies, resulting in 16,048 protein-coding genes. This resource provides a foundation for genomewide analyses of inbreeding, demographic history, and adaptive variation, and will support evidencebased in situ and ex situ conservation strategies for this threatened species.

The fungal crop pathogen Blumeria hordei, causal agent of powdery mildew on barley, presents life-history and epidemiological characteristics, as well as and selective pressures due to modern agriculture leading to expected sweepstakes reproduction, that is highly skewed offspring distributions. Using genome-wide polymorphism data and population genomics inferences, we aim to 1) infer the past demographic history and the strength of sweepstakes reproduction in B. hordei, and 2) quantify the contributions of these selective and neutral processes in the genome. An new inference method based on Neural Posterior Estimation and diversity and linkage disequilibrium statistics was developed and tested on simulated and B. hordei genomic data. We confirm that B. hordei exhibits a moderate sweepstakes reproduction (-parameter of 1.6). We highlight that the Site Frequency Spectrum (SFS) appears sensitive to the joint occurrence of sweepstakes and recent demographic changes, which may caution on the reliability of the SFS to infer sweepstakes reproduction. We then scan the genome for selective sweeps, adjusting the significance thresholds of the methods for demographic history and sweepstakes reproduction, thereby yielding a counterintuitive result. When conditioning the significance threshold for sweep detection on simulations under sweepstakes and demography, a very large number of putatively selected regions is found (11.6% of the genome). We suggest that sweepstakes reproduction in B. hordei is due to 1) neutrality (clonal/sexual phases and Boom-and-Bust cycles) generating a genome-wide level of background noise in the coalescent genealogies, and 2) selective sweepstakes due to pervasive positive selection. Our findings have important implications for both population genomic methodology and our understanding of pathogen evolution.